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^ . Abstract: We propose as a generalization of an idea of Ruelle to describe turbulent fluid flow 

a chaotic hypothesis for reversible dissipative many particle systems in nonequilibrium sta- 
00 ' tionary states in general. This implies an extension of the zeroth law of thermodynamics to 

non equilibrium states and it leads to the identification of a unique distribution fx describing 
the asymptotic properties of the time evolution of the system for initial data randomly chosen 
^ . with respect to a uniform distribution on phase space. For conservative systems in thermal 

lO \ equilibrium the chaotic hypothesis implies the ergodic hypothesis. We outline a procedure 

to obtain the distribution fi: it leads to a new unifying point of view for the phase space 
behavior of dissipative and conservative systems. The chaotic hypothesis is confirmed in a 
' non trivial, parameter-free, way by a recent computer experiment on the entropy production 

, fluctuations in a shearing fluid far from equilibrium. Similar applications to other models 

' are proposed, in particular to a model for the Kolmogorov-Obuchov theory for turbulent flow. 

>>■ 

O ! ^1 Introduction. 

a 

' In a previous paper [CG] we proposed the use of RueUe's idea (discussed in §2) to obtain 

the probabihty distribution for the statistics of turbulent flows in hydrodynamics, as a basis 
for the study of many particle statistical mechanical systems in nonequilibrium stationary 
^ ' states in general. We did so, by providing a concrete procedure of how to obtain the 

necessary probability distribution, now called the Sinai-Ruelle-Bowen (SRB) distribution, 
to compute the statistical properties of the above mentioned systems. The applicability of 
such a distribution has, so far, only been proved with mathematical rigor for very idealized 
systems, such as Anosov or Axiom A systems, and it would be impossible at present to give 
extensions of the proofs for the many particle systems of interest here. Therefore we proposed 
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to use Ruelle's idea as a heuristic principle to obtain the statistical properties of such systems, 
at least when they are very large, i.e. in the thermodynamic limit. This implied that we 
made a "chaotic hypothesis" that the many particle systems in statistical mechanics are 
essentially chaotic in the sense of Anosov, i.e. they behave as if they were Anosov systems 
as far as their properties of physical interest arc concerned. In other words we use the SRB 
distribution obtained from the strong assumption of chaoticity in the Anosov sense in a 
heuristic way to compute statistical mechanical properties of our system and assume that 
the corrections due to the possible non-validity of the strong chaoticity assumption become 
negligible for large systems. 

The position that we take here is very similar to that usually taken with respect to the 
so-called ergodic hypothesis, which has been proven only for very special few particle sys- 
tems. Yet, when used as a principle, it has led to all known results of statistical mechanics, 
beginning with its connection with thermodynamics. It would seem therefore inappropriate, 
in fact very unfortunate, if the application of the ergodic hypothesis, would have had to 
wait till it had been proved valid for the many particle statistical mechanical systems in 
thermal equilibrium of physical interest. Very recently a version of what we shall call the 
chaotic hypothesis, see §2, has been rigorously proved for a single particle system held in a 
non equilibrium stationary state and a number of detailed consequences have been derived, 
which agree with experiment [CELS] . 

Here we will give a number of possible many particle systems to which the chaotic hy- 
pothesis and the ensuing SRB distribution could immediately be applied. So far, only one 
of those systems: a shearing thermostatted fluid far from equilibrium (see §3, model 2) has 
been investigated experimentally, viz. the statistics of the fluctuations of the pressure ten- 
sor - or equivalently of the entropy production rate - in this system have been determined 
numerically and found to be in very good agreement with what one obtains by applying 
the chaotic hypothesis. Although corresponding experiments have not been done as yet for 
the other systems we mention, they should provide further checks on the validity of Ruelle's 
ideas and the chaotic hypothesis as proposed here. 

We want to emphasize that the application of the chaotic hypothesis is not restricted 
to stationary states near equilibrium, i.e. to the linear regime of small deviations from 
thermal equilibrium, as the above mentioned example of a shearing flow shows. The precise 
limitations of its applicability are unknown, however. 

The way we will present the construction of the SRB distribution from the chaotic hy- 
pothesis can also be applied to the theory of equilibrium states. It leads then to a new 
picture of the behavior in phase space of both equilibrium and stationary nonequilibrium 
systems, which reveals a much closer analogy in their phase space behavior than considered 
up till now. Thus a unification of the conservative behavior in equilibrium states and of the 
dissipative behavior in non equilibrium stationary states emerges. 

In §2 we describe some general properties that can help visualizing the general phenomenol- 
ogy of the non equilibrium systems that we consider: the discussion leads then to a formal 
definition of Ruelle's idea and to the precise formulation of the chaotic hypothesis. In §3 
we give a variety of examples of nonequilibrium systems to which the chaotic hypothesis 
can be straightforwardly applied. In §4 we discuss from a somewhat unusual viewpoint the 
heuristic ideas behind the hypothesis; this leads, in §5 and §6, to an outline and reinterpre- 
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tation of the classical, [S2,Bo,Rl], construction of the appropriate SRB distribution for this 
system, using Markov partitions. In §7 wc briefly summarize the only concrete application 
so far available, viz. that of a shearing fluid, and we discuss our main result, the fluctuation 
theorem of §7 (which gives a theoretical interpretation of the experiment). In §8 we give a 
discussion and outlook. 

%2 The SRB picture. 

For a convenient discussion of the SRB picture of nonequilibrium stationary states it is 
important to discuss the time evolution in discrete time, rather than in continuous time. 
This will be obtained by observing the motion when some timing event happens (this is 
usually done by describing the motion through a Poincare section). Therefore we fix a 
timing event and envisage performing our observations at every time the event happens. 
This will have the effect of reducing by one unit the initial phase space dimension. 
The choice of the timing event is essentially arbitrary: for many particle systems a reasonable 
choice could be the event in which the pair of closest particles (among the N we have) is at 
a distance r, coming from larger distances. Wc call such event a "collision" and wc use it 
as our timing event. To avoid trivialities r has to be chosen small compared to the average 
interparticle distance, but not too small {i.e. larger than the "core" of the interaction). In 
the case of a continuous fluid flow a natural timing event, actually used in many numerical 
experiments starting with [Lo], is the event in which a given coordinate of the velocity fleld 
passes through a prefixed value, or assumes a locally maximum value. To uniformize the 
notations we shall also call such an event a "collision" . 

The dynamical systems we consider will be defined now by the phase space C of the "colli- 
sions" , with dimension 2D and the time evolution, which will be a map 5' : C — > C defined by 
Sx = if St is the continuous time evolution operator solving the equations of motion 

in the full phase space J-', which in our cases will coincide with a constant "energy" surface 
or it will be a manifold in it. Here t{x) is the time interval between the collision x G C and 
the next one. 

We shall make a statistical study (like is done in equilibrium) : this means that we shall be 
interested in the properties of the time evolutions of the motions that can be seen by extract- 
ing the initial data at random with the Liouville distribution on J^. Since our analysis will be 
performed on C rather than on we shall need the corresponding probability distribution on 
C. The Liouville distribution fj,c = const 5{H{p, q) — E)dpdq (or /X£ = const |7^|^ — E) 
in the case of the fluid motion models that we consider in §3, where the variables 7^ are the 
Fourier components of the velocity field) on the full phase space !F (energy surface) naturally 
generates a probability distribution /xq on C: if is a set on C we simply set lJio{E) equal to 
the Liouville measure of the tube of trajectory segments in T starting at E and ending at 
the next collision, when evolved with the motion corresponding to no external forcing fields. 
We shall still call /xq the "Liouville distribution" (on C). 

The first point of our analysis is a generalization of the zeroth law of thermodynamics to 
nonequilibrium stationary states. As expressed by Uhlenbeck and Ford, [UF], the zeroth 
law of thermodynamics states that a closed conservative mechanical system consisting of a 
very large number of particles will, when initially not in equilibrium, approach equilibrium, 
where all macroscopic variables have reached stationary values. By (asymptotic) equilibrium 
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one means here that the time averages have reached the value that can be computed by a 
probability distribution on the energy surface. This law can be extended to nonequilibrium 
systems as: 

Extended zero-th law: A dynamical system {C,S) describing a many particle system (or a 
continuum such as a fluid) describes motions that admit a statistics n in the sense that, 
given any (piecewise smooth), m,acroscopic observable F defined on the points x of the phase 
space C, the time average of F exists for all jjLQ-randomly-chosen initial data x and is given 



where n is a S -invariant probability distribution on C. 

In this form we suppose that it holds for all our models. The notation ii{dx) rather 
than r{x)dx expresses the possible fractal nature of the support of the distribution /z, and 
implies that the the probability to find the dynamical system in the infinitesimal volume dx 
around x may not be proportional to dx, so that it cannot be written as r(x)dx with r{x) 
a probability density and dx the volume measure on phase space. 

The main point of this paper is to use an idea of Ruelle's as a guiding principle to describe 
nonequilibrium stationary states in general. That is, we propose that for such systems there 
exists a distribution (usually called the SRB distribution) describing the asymptotic statistics 
of motions with initial data randomly chosen with respect to a uniform distribution on phase 
space (the Liouville distribution). For this to be realistically implemented wc assume that 
macroscopic systems, consisting of very many particles, behave as transitive Anosov systems, 
i.e. are "chaotic" in the sense that each point x in phase space admits an unstable and a 
stable manifold W^,W^ which depend continuously on x, are dense in the phase space C, 
and on which the expansion and contraction rates are everywhere separated by a finite gap 
from (hence no zero Lyapunov exponents occur). ^ 

We propose therefore the following chaotic hypothesis, which in [CG] we called Ruelle's 
principle, as a generalization of Ruelle's idea: 

Chaotic hypothesis: A reversible many particle system in a stationary state can be regarded 
as a transitive Anosov system for the purpose of computing the macroscopic properties of 
the system. 

We intend to show that this hypothesis, although general, leads to concrete verifiable 

^ For convenience we formally recall that an Anosov system (C, S) is a smooth dynamical system such that 
every point a; G C possesses stable and unstable manifolds , which depend continuously on x and on 
which S",S~" respectively contract infinitesimal vectors by a factor bounded by Ce~^", for n > 0, and 
likewise for n < they expand by a factor bounded by C~^e~^". The constant A is therefore such that all 
Lyapunov exponents verify |Aj| > A and, hence, A can be called a bound on the Lyapunov spectrum gap. 
Note that the continuity of the W^, in x implies the transversality of the two manifolds, which therefore 
form everywhere an angle '&{x) bounded away from and tt. An Anosov system is transitive if W^, are 
dense in C for all x. 



by: 
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conscqucmccis and may be, in this respect, similar to the ergodic hypothesis for equiUbrium 
states but, unhke the ergodic hypothesis, admits an extension to nonequiUbrium stationary 
states. 

One could weaken our form of the chaotic hypothesis by replacing "Anosov system" with 

"Axiom A system" , and refer to the general theory of such systems developed in [Bo] . [R2] 
(rather than relying on the work of Sinai on Anosov systems, [S2]); one could even attempt 
to weaken it further by trying to make use of the general theory of Pesin of non smooth 
hyperbolic system, [P]. However, we shall not dwell on such somewhat obvious extensions 
of our ideas, as they do not seem relevant at present. 

Examples of model systems in nonequilibrium stationary states to which the chaotic hy- 
pothesis is applicable will be given in §3. Numerical evidence leads us to believe that they 
seem to share a number of properties which we believe to hold also for more general physical 
systems and which we now summarize. Not all of them are necessary for the applications we 
shall discuss: however, they are very helpful for building an intuitive, model independent, 
picture of the phenomena that we attempt to study. When discussing the models from a 
technical viewpoint we shall mention which properties have been experimentally checked and 
which have not (yet) been checked: our applications will only require properties (A,B,C) 
below. 

(A) Dissipation: the phase space volume undergoes a contraction at a rate, on the average, 
equal to D{a{x))_^_ where 2D is the phase space C dimension and a{x) is a model dependent 
"rate" per degree of freedom. The average here is a time average from time zero to plus 
infinity and the rate is a generalization of the usual entropy production rate, (see §3 for 
motivation of this remark). 

We say that a system is dissipative if the contraction rate per degree of freedom, {a-)_^_, is 
positive. We shall assume that the models that we consider here in nonequilibrium situa- 
tions are all dissipative. The instantaneous contraction rate ^{x) is, however, a fluctuating 
quantity and we note that when we consider in this paper entropy production rates and 
their fluctuations we identify them, mathematically, with phase space contraction rates and 
their fluctuations, respectively. 

(B) Reversibility: there is an isometry, i.e. a metric preserving map, i in phase space, which 
is a map i : x ix such that if t — > x{t) is a solution, then i{x{—t)) is also a solution, and 
furthermore is the identity. 

(C) Chaoticity: the above chaotic hypothesis holds and we can treat the system {C,S) as a 
transitive Anosov systems. 

We realize that (C) cannot hold strictly in finite systems, even in the case of smooth 
interaction potentials (in the presence of hard cores the Anosov property, which requires 
smoothness of the dynamics as a prerequisite, is in fact obviously false). What we mean 
here is that we assume that the system behaves as if it was a transitive Anosov system and 
that the errors made become negligible (even when there are hard core collisions) at least 
in the large system limit. 

In support of (A,B,C) the following two properties (D) and (E) also are relevant and 
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appear to hold at least for some of the models that we shall treat: 



(D) Pairing of Lyapunov exponents: half of the 2D Lyapunov exponents are > and half 
are < 0. If they are ordered so that < Aj^ < < . . . < and > AJ~ > A^ > . . . > A^ 
so that Amax = A^ the following pairing rule holds: 



which has been proved for some special cases where the system is non reversible {i.e. a{x) 
is constant) and the pairs do not necessarily consist of exponents with opposite sign; it 
has been found numerically, in the case of model 2 and related models in the form (2.2) in 
[ECMl], [SEM] (where it was formulated in the above form). 

On the basis of what is presently known, one can conjecture that even if the pairing rule 
does not hold in the above form it could still hold in the form of an inequality: — (<t)_,_ < 
At + A J < {weak pairing rule). 

(E) Smoothness of the Lyapunov spectrum: the Lyapunov exponents become for large A'' a 
smooth function of their index. This means that, with the labeling of the exponents as in 
(D) above, if one draws a graph of x = — > A+ = fN{%), then in the "thermodynamic 
limit" (AT — > 00 with constant density for particle systems; in the case of fluid systems 
the role of N will be taken by the Reynolds number) fN{x) — > foo{x) where foo{x) is a 
smooth increasing function oi x G [0, 1] . 

Evidence for the generality of such property comes from [LPR], and [ECMl], [SEM], and 
quite likely it holds for all the models we consider in §3. 

We make the following remarks on the properties (A) to (E). 

1) First we note that the irreversible entropy production {a)_^_ in (A) results in a phase 
space volume contraction. This implies in turn that the attractor that we denote (which 
by property (C) is just the full phase space C) will contain an invariant set A of zero Liouville 
measure and dimension equal to the fractal dimension of the motions (and strictly less than 
that of the phase space, see below) but of probability 1 with respect to the statistics of the 
motions generated by the dynamics S from initial data chosen randomly with respect to the 
Liouville distribution /iq. 

It is convenient, therefore, to distinguish between the attractor Aq as it is usually defined in 
the literature (which is a closed set for virtually all adopted definitions) and our sets A. The 
latter are not uniquely defined, but they are in an obvious sense more intrinsically related 
to the motions. It can very well be that Aq is smooth and even coincides with the full phase 
space, as is the case when (C) holds, while A is much smaller (and is a fractal). Thus in 
this paper, unlike in most established conventions, we shall call A the attractor. it will not 
matter which particular A one considers. 

We adopt, as definition of the fractal dimension of the motions {i.e. of A), the Kaplan 
Yorke definition (also called the Lyapunov dimension, [ER]) The latter is, probably, [ER], 
quite generally equal to the Hausdorff dimension of those sets A which have the smallest 





(2.2) 
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HausdorflF dimension and which are visited with frequency 1 by almost all, with respect to 

the Liouvillc distribution /io, motions [ER] p. 641. 

2) The above properties (D,E) imply that the attractor A for the motions with a given 
energy is a fractal set with a dimension close to the full dimension 2D: the fractal dimension 
will be, in fact, of the order of 2D — 0{{a) ^X^\^)D as immediately follows if one adopts, 
as above, the Kaplan Yorke definition of fractal dimension. Note that (E) and the above 
weaker pairing rule are sufficient for this conclusion. Systems for which smoothness and the 
(weak) pairing rule hold do show dimension reduction, i.e. the attractor in phase space has 
a dimension which is macroscopically different from that of the phase space itself. 

3) Reversibility implies an important property of the attractor A: if A is an attractor 
for the forward motion then A_ = iA is an attractor for the backwards motion and, more 
generally, the statistical properties as t ^ ±oo, of the motions generated by initial data 
randomly chosen with respect to the Liouville distribution no are trivially related. 

4) The basic properties for the validity of our results for nonequilibrium stationary states 
are chaoticity (C) and reversibility (B). If (C) holds the existence of the SRB distribution, 
i.e. of a probability distribution describing the asymptotic statistics of the motions of a 

system evolving with a dynamics S, whose initial data are chosen randomly with respect to 
the Liouville distribution hq in the phase space C, can be proved as a theorem: 

Theorem: if a system {C, S) is a transitive Anosov system then it admits a SRB distribu- 
tion, [S2J. 

In conservative systems in equilibrium, verifying (C), the distribution /x is the same as the 
Liouville distribution /xq itself (which is invariant by the Liouville theorem): see [S1,S2],[AA]. 
Hence (C) implies the ergodic hypothesis in this case and the attractor A can be taken to 

be the full phase space C. 

In this paper we are interested in dissipative systems satisfying (A) where new phenomena 
occur and ^ is not the Liouville distribution fiQ. For systems verifying (A),(B),(C) we prove 
a fluctuation theorem, see §7, which is our main technical result. 

§5 Models. 

We now list a number of models to which our theory can conceivably be applied. All 
these models contain thermostat mechanisms in order to enable the systems to reach a non 
equilibrium stationary state in the presence of an imposed external field. Model 1 is a model 
related to electrical conductivity, models 2 and 3 are related to shear flow, model 4 to heat 
conduction and model 5 to a fluid mechanics model for turbulent flow. 

We distinguish, as in §2, between the phase space over which the system evolves ac- 
cording to the equations of motion and the collision phase space C consisting of the timing 

events ("Poincare section of JF"). 

The details of the models described here will not be used in the following since our main 
point is the generality of the derivation of a fluctuation formula from the chaotic hypothesis 
and its (ensuing) model independence. However, we include them for concreteness and 
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reference. 



Model 1: a gas of A'' identical particles with mass m, interacting via a stable short range 
spherically symmetric pair potential ip and with an external potential (p'^ ^ 0, enclosed in 
a box [—\L, and subject to periodic boundary conditions and a horizontal constant 
external field Ei {i is a unit vector in the a;-direction). The external potential will be just 
a hard core interaction which excludes access to a number of obstacles (hard disks, to fix 
the ideas) so situated that every trajectory must suffer collisions with them. The system 
is in contact with a "thermostat" adding (or subtracting) energy so that the total internal 
energy stays rigorously constant. The equations of motion are: 

1 

ij = ^Pf E,=F,+Et-a{p)p^ 

J ^ — / _j —J — ( _j —J 

with j = 1,...,N; a{p) = Ei - J2jPj/ {J2j P^) ^^'^ flj is the force acting on particle j. 
The a-term incorporates the coupling to a "gaussian thermostat" and follows from Gauss' 
"principle of least constraint" . The constraint here is the constancy of the internal energy: 

^(£'£) = E + E ^(«. - + E ^'(li) (3-2) 

i=l i<j i 

a typical nonholonomic constraint; it follows then from Gauss' principle that the force cor- 
responding to the constraint is proportional to the gradient with respect to p^ of H. This 
model has been studied in great detail, in [GELS], in the case = 1 and <f = and ip^ a 
hard core potential as above, making it a Lorentz model for electrical conductivity if E is 
an electric field; a similar model has been investigated numerically in [BEC]. It is part of a 
wide class of models, together with the following models 2,3, whose interest for the theory 
of non equilibrium stationary states was pointed out in [HHP], [PH], where one can find 
the first studies performed in the context in which we are interested. The dimension of the 
phase space T of this system is do = 4:N — 1 and that of C is do — 1 = 2D with D = 2N — 1. 
The phase space "contraction" rate, i.e. the divergence of the right hand side of (3.1), is 
Da{x) = Da{x), which can be written in the form: 

Dai.)=Dai.)=D^^ (3.3) 

where e{x) is the work done on the system per unit time by the external field and kT{x) is 

h m which, if k is Boltzmann's constant, defines a kind of kinetic temperature: hence 
the name of entropy production rate per (kinetic) degree of freedom that will be occasionally 
be given to a{x).^ Note that cr{x) does not have a definite sign. 



^ If p is the average of i ■ p . , then i {p^, ) = + mkT and we see that T{x) cannot be identified with 
the temperature unless we neglect p^ compared to mkT, i.e. we identify the peculiar momentum, needed for 
the proper definition of the temperature with the ordinary momentum (which should not be done at large 
E where one cannot identify p^) with mkT, with T being the usual temperature). 
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It has been proved, [CELS], that for N = 1 and smaU E > the average (cr)^ is positive, 
i.e. the system is dissipative in the sense of §2. There seems to be no reason to think that 
(a)_|_ is not positive. For the above model with N > 1 no experiments are available yet on 
the pairing rule or the Lyapunov spectrum smoothness. Nevertheless one can present an 
argument for the validity of the pairing rule, which may sound convincing but that we have 
been unable to substantiate mathematically."^ In this case the time reversal map i is just 
the usual i : {q,p) (q, —p). 

Model 2: a shear flow in a two dimensional container [—^L, where the particles evolve on 
a moving background running with velocity vyy in the x direction proportional to the height 

y in the y direction, is measured with respect to that of the center of mass of the particles. 
The background exercises a drag on the j-th particle located at height yj proportional to its 
peculiar velocity: q . — i^yj with respect to the background. The introduction of yj instead 
of the usual yj is due to the boundary conditions we choose (see below), which do not keep 
the height of the center of mass of the particles fixed. For large N the difference between 
yj and jjj will become negligible (see comment 6) in §8). Similarly for large TV (and large 
L with n = NL~^ fixed), the forcing 7 should be the shear rate in the fluid, i.e. 7 = 
where Ux{y) = i"fy is the average (local) velocity of the particles in the fluid at height y, 
as is indeed found in the computer experiments for this system, [ECM2] . If F_j is as in the 
second equation of (3.1) and (fi^ = and q. = {xj,yj), the equations of motion are, using 
Gauss' principle of least constraint to keep the internal energy fixed: 

i2. = ^^^"(^,'i7y,)a ^^^^ 

with j = 1,...,N. Here p. = m{q . — v^yj) is the peculiar momentum relative to the 
background flow; F_j ^ is the force on particle j due to particle i and a is again defined so 
that (3.2) is a constant of the motion; finally 7 plays here the role of a forcing field as E 
did in model 1. One imposes periodic boundary conditions on the horizontal direction; on 
the vertical direction a natural boundary condition is perfect reflection against the walls at 
y = ±i_L. This model has been extensively studied, numerically, in [ECMl], [ECM2] with 
somewhat different boundary conditions. 

One can suppose that the total horizontal component of the peculiar momentum and 
the horizontal position of the center of mass denoted, respectively, P^tX^., are 0: this 
is consistent with the equations of motion. We shall refer to Px and as conserved 
quantities: but one should bear in mind that they are such in an "improper" way because 



Suppose that the equation (3.1) is modified into the same equation with a(p) constant. Then [Dr],[ECMl] 
prove that the 4^ Lyapunov exponents can be paired so that the sum of the corresponding pairs is just 
—a. In model 1 a is not constant: however it has an average value (a) which is constant on the attractor, 
with probability 1 with respect to the choice of the initial data (with distribution ^0): therefore we may 
hope that "things go as if" a was constant; hence the Lyapunov exponents should be so paired that their 
sum is —{a). This is not yet the above full pairing rule because there we assert in addition that half of the 
exponents axe positive and half are negative: and this only "follows" if reversibility (B) is also used. 
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they are conserved only if their initial values are 0: P^; in fact relaxes to with a Lyapunov 
exponent which in general is not zero (and equals the time average (a)_,_ of a). 

If it is assumed that Px = = Q and if one recalls that also is a constant of the motion 
and one imposes a priori its value, then the dimension of the phase space ^ is c?o = 4A'' — 3, 
which wc write c?o = 2D + 1 for imiformity with the notation in model 1, so that D = 27V — 2. 
The phase space contraction rate, i.e. the divergence of the r.h.s. of the equation of motion 
regarded as first order equations for p, q is Da{x) with: 

a{x) = a{x) + 7§g^ = a{x) + ^0{N-^) = (3.5) 

where T{x) can actually be interpreted as a kinetic temperature, so that (7{x) can also be 
called the entropy production rate. 

Numerical experiments with N up to 864 show that ((t)_^ > 0, [ECMl], ECM2], [SEM]. 
Such papers also provide (strong) evidence for the pairing rule and some (weak) evidence 
for the smoothness. The time reversal map in this case is not the usual velocity reversal but 

1 : {x,y,Px.Py) {x, -y, -Px,Py)- 

Model 3: this is a model for a shear flow, produced by boundary forces, in contrast to model 

2 where the shear is produced by body shear forces. 

The flow proceeds in a two-dimensional container [—■§,■§] and the equations of motion 
are simply: 

q. = —p., P.=F^ (3.6) 

supplemented by periodic boundary conditions on the horizontal direction and shear gener- 
ating boundary conditions in the vertical direction: 

= m (3.7) 

where w is the collision angle formed by the incoming velocity with the x-axis, counted coun- 
terclockwise for collisions at y = L and clockwise for collisions at = 0; u' is the outgoing 
velocity angle formed with the a;-axis, counted clockwise at y = L and counterclockwise at 
y = 0. 

With the above angular conventions, ui' = u represents the ordinary elastic collision. We 
shall consider a "shearing collision ndc" uj' = f{u}), uj' < oj, where / is a reversible collision 
rule. Reversibility here has the literal meaning: the collision obtained by reversing the 
particle velocity after a given collision {i.e. the incoming collision with an angle tt — /(w)) 
produces afterwards the reverse of the original collision, {i.e. it — u). This is: 

7r-u; = /(7r-/(a;)), fico) < uj (3.8) 

where the first condition is the reversibility condition and the second is the shearing con- 
dition. The equations (3.8) can be solved by simply thinking that the graph of f{ui) as a 
curve u — > {u),f{u))) e [0,7r]^ is a concave arc connecting the point (0,0) with the point 
(tt, tt), symmetric by reflection around the secondary diagonal of [0,7r]^. Furthermore one 
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imposes that the colhsions preserve also the horizontal total momentum and the total en- 
ergy. Since the horizontal momentum of the colliding particle (say particle 1) changes (by 
\p^\{cosu)' — cosw)) one can impose the two conservation laws by gaussian minimal con- 
straints, i.e. by requiring that the variation of the other momenta is: 

p;. - (1 + 'd)p^ + Pi (3.9) 

for suitable multipliers d and /3, which after a brief calculation leads to explicit expressions 
for /3, 1?, with /J = 0{N~^) while •& = 0{N~^), so that the relations (3.9) generate corrections 
which, while enforcing the constraints, can be regarded as negligible for large N (see comment 
6) in §8). 

The symmetry of / also guarantees that the collision rule (3.8) corresponds to a gaussian 
constraint (forcing \ to stay constant in the collision, while changing the direction with 
respect to the elastic collision), [CL]. 

For concreteness one can take, following [CL], / such that the above arc is an arc of a 
circle centered on the secondary diagonal and passing through the indicated points. The 
circle curvature will be a measure of the shear strength. The dimension of the phase space 
of this system is do = 4:N — 3, if we fix the energy, the horizontal total momentum and 
the horizontal position of the center of mass. We write do = 2D -|- 1 as in the previous 
models so that the dimension of the collision space C is 2D with D = 2N — 2. We also 
suppose, naturally, that the collisions with the walls are among timing events. Then at 
every collision there is a reduction of phase space volume = ^"^[^^ f'i'^): [CL]- If 

we define n{x) = 1 when a; is a collision with the wall and n{x) = otherwise, the phase 
space contraction can be conveniently written as e^^"^*^^^""'^^ with t{x) equal to the time 
elapsing between the collision at x and the next at Sx and with the entropy production rate 
a{x) defined by: 

Daix) = ^ log(l - ^^/'(.))- (3.10) 

This model has been studied in detail in [CL].^ 

There is numerical evidence that if fiuj) ^ uj then (cr)^ > 0. The pairing and smoothness 
properties have not yet been studied. The time reversal operation is the "usual" one (see 
model 1). 

Models 1 and 2 arc easier to interpret as dynamical systems than this model, but they are 
physically somewhat artificial, in that the gaussian thermostat is a rather unconventional 

^ For this purpose one has to think, [CL], that the coUision rule (3.8) is realized as a limit of a gaussian 
constraint rule acting on a tiny corridor of width S expanding vertically the container at the top and bottom, 
where the particles can enter from the inside but are subject to a horizontal field ±£5 (the sign depending 
on whether the particles are in the upper or lower corridor): the field will produce a bias in the scattering 
angle such that the particle will come out with an angle different from the incoming angle. The constraint is 
that the kinetic energy of the particles inside the corridors does not change; the particles colliding with the 
external corridors walls are just perfectly reflected. In the limit as 5 — > and — ► 00 (at suitable rates) a 
reflection rule like (3.8) is realized with a special /. By letting depend on the distance to the corridor 
boundaries essentially any / can be realized, in the limit 5 — ► 0. 

^ It is not difficult to see, by thinking of the constraint as a (limit of) gaussian constraints as above, to see 
that also in this model Da{x) is an "entropy production rate" . 
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model for a thermostat and the shear force is a body force rather than the usual boundary 
force. In this respect model 3 is better as its unphysical reflection laws are a boundary effect 
which only produces the net effect of generating a shear on the system. 

Model 4) a model for heat conduction considered in [HHP], [PH]. In a box, [—L — 
H,L + H] X [— "I], N particles move interacting via a short range pair potential; the 
boundary conditions are perfect reflection horizontally and periodic vertically; the particles 
are subject to the nonholonomic constraint that the total kinetic energy in the left part of 
the box [—L - H, -H] x [-■§, f] and in the right part of the box [H,L + H] x [-^, ^] have 
constant values, denoted, respectively by 2l+2h ^^'^- ^^'^ 2l+2h ^-^+' — obey: 

N p'^ N T 

i=i 

Here x± E^re the characteristic functions of the left and right parts of the box which have to 
be interpreted as plates of thickness L at temperatures T+ and T_ respectively; q = {x,y). 
In addition wc impose that the total energy (given by (3.2) with <fi'' = 0) of the gas = H 
is exactly conserved. 

The constraints are implemented using Gauss' principle, i.e. by a force proportional 
to the gradients with respect to the p^.'s of $± and $o (which are simply x±(a;j)^, 
respectively), leading to the equations of motion: 

gj = ^Pj^ Ej = £j ~ oi+X+ixj)Pj - a-X-{xj)Pj - aop^ (3.12) 

with q;±, ao are defined so that <1>±, $o are exact constants of motion. The values of a±, ao 
can be easily computed; their general expression will not be needed here. For the purpose 
of illustrating once more that the resulting forces will lead to a reversible dynamics wc give 
their expression in the simple case in which only $± are imposed: in this case the values of 
a± are relatively simple and ao is not present. One finds: 

"± = ^.(^.\2^2 (3.13) 



Going back to (3.12) we note that, with the three mentioned constraints, model 4 should 
be a quite realistic model for heat conduction. The dimension of the phase space J- is 
do = 4:N — 3, which again we write as do = 2£> + 1, with D = 2N — 2, so that the dimension 
of the collision space C is again 2D. The phase space contraction rate is in this case: 

Daix) = a+{x)2N+ + a_(x)2iV_ + ao{x)2N + O(A^-i) (3.14) 

if N± denote the number of particles in the right and left "plates". Eq. (3.14) could also 
be interpreted as in the previous models as an entropy production rate. Some numerical 
evidence that (cr)^ > if > T_ can be found in [PH],[HHP]. No evidence for pairing or 
smoothness rules seems available. The time reversal map is the "usual" one, see model 1. 

The above model equations can be made smoother by replacing % by a smoothed version 
of the characteristic functions of the plates; say by functions which are = 1 except within 
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a distance of the order of the interaction range from the inner boundaries of the plates: in 
this region x± decrease gently from 1 to zero. 

Model 5): this is a model related to turbulent flow, obtained from the Navier Stokes 
(NS) equations. We consider the NS equations in a box [— -j]^, with periodic boundary 
conditions and for an incompressible fluid. If the velocity field is written in a Fourier series 
as: 

«fe)=E^--Ifc (3.15) 

with 7^ complex vectors with 7^ = 7_^ (reality of the velocity field) and 7^ ± fc (incom- 
pressibility) then the NS equations become: 

ik = -' E (lfe/^2)%7,^+i?4-^^fc^, (3.16) 

Hk is the orthogonal projection over the plane orthogonal to k: v is the kinematic viscosity 
and R.f ^ is the forcing (of course orthogonal to k) which will be taken to be non zero only for 

a few components with small k. Since k = with n integer this means that the force acts 
only on the high length scale components. For simplicity we may think that the forcing has 
only two non vanishing components -R/^o 1 -R/j,o j corresponding to two linearly independent 

wave numbers ki.k2-^ The number R therefore determines the forcing strength and will 
be identified with the Reynolds number (we keep the container size L and the viscosity u 
fixed). We take 7^ = since it is the conserved center of mass velocity. 

In order to obtain equations in the framework of this paper from the phenomenological 
theory of Kolmogorov-Obuchov, [LL], we shall assume that the above equations can be 
replaced by the following simpler ones: 

ik = -' E_ dfc, • ^2) nfc7fe, + 4 

Here, if ka = the wave vector kfj is the Kolmogorov momentum scale k^ = k^R^/'^, ([LL] 
p. 122, (32.6)), so that if Nr is the number of wave vectors ("modes") k such that when 
fco < |fc| < fcfl + v~^/'^ then A''^ k. (^)^ ~ i?^/^; then the phase space C has dimension 
2Nr -2 = 2D with D « R^/'^, while has dimension do = 2Nr -iJ 

This means that the equations for the amplitudes 7^ corresponding to fc's in the inertial 
range, fcg < |fc| < kn, are "governed" by the reversible Euler equations. In the viscous range, 

^ The simpler case of only one non zero component can be trivial {e.g. if the forcing acts on the smallest k, 
\k\ = ko) and is therefore discarded here in favor of the next to the simplest, [Ma]. 
Taking into account the reality and incompressibility conditions, forces the 7^ to have only two linearly 

independent components. 



|fc| < fcj 



kR < \k\ <kR + 



(3.17) 
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|fc| > ku the dissipation phenomena wiU be idcahzcd by saying that the equations are simply 
such that only the modes k with ka < \k\ < ka + v~^l'^ have a non zero amplitude and 
evolve in such a way as to keep the total energy constant. This means that the parameter 
a is an effective thermostat (or viscosity), which has to be chosen so that the total energy 
is constant, i.e. so that ^ |7j_p = 0: 

Here e{x) and DykT{x) are simply the numerator and denominator of the fraction defining 
a{x), if 2Nyii is the number of modes in the viscous range and one defines 2D„ = 2Nyji. 

The Kolmogorov length k^^ is introduced here phenomenologically and we do not attempt 
at a fundamental derivation of (3.17), (3.18). Therefore (3.17) has to be regarded as a 
phenomenological equation. 

Note that a is proportional to the work s{x) per unit time and per viscous degree of freedom 
performed on the system, which is dissipated into heat, in order to keep the total energy 
constant: the proportionality constant is 2DykT{x) with kT{x) = X]|fe|<fcji It^P (which, 
however, is not a constant of motion for (3. 17), (3. 18) because of the imposed constraint that 
X]fe It^^I^ rather than '}2\k\>kR I'^fel^ constant). The phase space contraction rate is in this 
case: ~ ~ 

D,a{x) = D,a{x) = D,^^^^ (3.19) 

Hence Dy{a)_^_ can be thought of as the average amount of energy dissipation per unit time 
by the flow divided by the kinetic energy contained in the viscous modes. The first quantity 
plays a major role in Kolmogorov's theory, see [LL] p. 119, and its average is usually called e, 
(see [LL], (31.1)). Since the kinetic energy contained in the viscous modes can be thought of 
as a kind of "temperature" we see that 2Dy{a)_^_, is proportional to the entropy "production 
rate". More appropriately we can say that, for R large, 2Dy{a)_^_ is proportional, once 
more, to the "energy dissipation rate" over a kinetic quantity equal to the average kinetic 
energy contained in the viscous modes if, for large R, the two quantities can be regarded as 
independent random variables. 

Note, however, that for the above model (3.17) (introduced, we believe, for the first time 
here) there is no evidence for (cr)^ > or for the pairing and smoothness rules. The time 
reversal map is simply i : {7^} {~7^}- 



The introduction of the SRB distribution. 

We now present a heuristic argument providing, in our opinion, a useful characterization 
of the SRB distribution: this point of view is important for the applications of §7. Our 
purpose is to look at it from a somewhat different perspective than usual and to show that 
it leads to a new interpretation of the ergodic hypothesis and to a unification of equilibrium 
and nonequilibrium statistical mechanics. 
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We deal with systems of N particles verifying the properties (A,B,C) of §2 and we observe 
their motions x S^x at discrete times in the collision space C of dimension 2D, see §2. It 
will be very useful to keep in mind the paradigm of hyperbolic systems: namely the Anosov 
map of the 2-dimensional torus T^: 



which plays a role analogous to that of harmonic oscillators in classical mechanics. This 
example is not only enlightening, but it is really the main source of intuition. Note that this 

is a reversible map if i is defined as i : {ipi, ip2) — > {ip2, ~'Pi), because iS = S~^i. 

Let O be a fixed point on the attractor A and let Wq be the unstable manifold Wq of O 
(dense on A, see (C), §2).^ The dimension of Wq is D, half that of the phase space C, see 
§2- 

For simplicity we shall suppose that O is a time reversal invariant fixed point O = iO; 
this assumption could be easily relaxed.^ 

The key idea on which we base our analysis is that the attractor A should be considered to 
consist of the smooth D-dimcnsional unstable manifold Wq of O (or of any other fixed point 
or periodic orbit in A with dense stable and unstable manifolds). Of course the manifold Wq 
can only fill A densely: we "lose" the accumulation points. But all the information needed 
to perform time averages should be already contained in Wq itself, as we are only interested 
in the averages of rather regular observables {e.g. picccwise smooth). We think that Wq 
coincides with A in the same sense in which the rationals can be regarded as coinciding 
with the reals in integration theory (which works only if one considers integrals of smooth 
functions) and can be used to compute numerically the integrals of smooth functions. In the 
same way statistical averages with the distribution (2.1) should be computable by simply 
approximating them with integrals over finite parts of Wq, like the sets ^-^A obtained by 
"blowing up" with a large time iterate a small connected surface element A of Wq 
containing O. 

In other words we want to regard the possible fractality of A as a rather irrelevant accident. 
We want to think of A as essentially identical to WqI the latter surface folds over and over 
again, being enclosed in the hounded phase space C. It therefore folds itself in C just as an 
uncut folio is folded into a book, thereby generating an almost three dimensional fractal set 
out of a two dimensional smooth manifold. But thinking of A as an unfolded manifold of 

Q 

There might be no such point fixed point O. However a periodic orbit starting at a point P and with 
period n would be a fixed point for 5" and we could get all the following conclusions by replax;ing S with 
5" , since the statistics of 5" and that of S coincide, when S is chaotic enough. Thus assuming the existence 
of a fixed point is not restrictive. 

^ It is not difficult to realize that in the models 1,2,3,4 there are always periodic orbits which are time 
reversal invariant, i.e. such that iO is also on the orbit; at least if one is willing to limit the particle density 

in some interval (whose size may depend on the range of the interaction). Also for model 5 it is very 
likely that periodic motions (unstable, of course) do exist. Note that since we are assuming (C), §2, it is 
automatically true that there are periodic orbits {i.e. chaotic systems always have a lot of periodic unstable 
orbits), [S2]. 




(4.1) 
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half the space dimension [D in our notation) leads to a change in the usual point of view, 
which regards ^ as a fractal set with dimension close to 2D. 

Introducing forcing and friction [i.e. passing from an equilibrium to a stationary non 
equilibrium problem) should then not be thought of as a real "discontinuity" : which would 
be the case if one took the viewpoint that one is passing from a nice smooth 2D dimensional 
attractor A = C to & nasty, strange, fractal A d C with dimension 2D — 0(q:A~J^^)Z), macro- 
scopically different from 2D (as implied by the pairing (D) together with the smoothness 
(E), §2). 

Rather it should be viewed as an insignificant deformation of the mistablc manifold Wq 
which will fold itself in C not exactly as in the conservative case, but leave a few holes 
between the "pages" to account for its global fractality. This is a change with respect to the 
conventional point of view for the case of conservative systems: these are no longer really 
different from the dissipativc ones. Their attractor has, in the new (unconventional) sense, 
exactly half the dimension of the full phase space (the dimension the conventional point of 
view attributed to them is that of the full phase space, i.e. twice as large). 

The main consequence of such a viewpoint, besides the mentioned unification of conserva- 
tive and dissipativc dynamics, is that it allows us to think of the attractor as "unfoldablc', 
with the consequence that our intuition about the motion on the attractor is greatly en- 
hanced. 

This unfolded attractor, imagined as a flat infinite surface, attracts exponentially fast 
nearby points: the approach to the attractor follows the stable manifolds associated with 
the attractor points, which can be thought of as needles sticking out of the attractor itself. 
The motion essentially consists, therefore, of an expanding [i.e. as unstable as possible) 
motion on the unstable manifold Wq. 

We can now easily understand the statistics /x, (2.1), on the attractor, i.e. the SRB 

statistics, as follows. 

Let [/ be a sphere with small radius h, centered at the fixed point O; and let us ask 
how to compute the time average of an observable F, if the initial data are chosen in U 
with uniform distribution, say, with a distribution absolutely continuous with respect to the 
Liouville distribution. 

Clearly the average of F over a large time T will be computable by looking at the image 
S'^U under S'^ for large T and by imagining S^U covered by the density into which the 
initial uniform density in U evolves in T time steps. If we call A the connected part of 

Wq n U this also means that we can regard the S"^ image, S'-^A, of the connected part of 
Wq nU as a good Enitc approximation A to our attractor. 

The set S'^U will be extremely thin and it will ''coaf the extremely large portion of Wq 
defined by 5^ A, [i.e. by our "good finite approximation" of the attractor), if we regard the 
attractor as unfolded. 

Let dx be a surface clement on Wq and let us regard S' as a map of Wq into itself. We 
shall call Au{x) the absolute value of the jacobian determinant duS (x) of S, as a map ofW^ 
into itself, at the point x. In this way A„(a;) will be the absolute value of the determinant of 
a matrix with a dimension equal to that of W^, i.e. D. Then we are interested in computing 
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the integral: 



/ pT{x)F[x)dx 



(4.2) 



where pT{x)dx is the amount of mass in the cyhnder with base dx, which is the image of the 
cyhnder in U with base S~'^dx and height equal to the height h of the initial "cloud of data" 
U . Denoting ds (resp. d^) the dimension of the stable (unstable) manifold of O (which in 
our case are dg = du = D), this means that pT{x)dx is proportional to h'^''\S~^ dx\, where 
\S~'^dx\ is the surface area of S~^dx. By the above definition of the local expansion rate 
A„(a;), one has then: 

-1 

PT{x)dx = const {S'"^ x)... A.~^ {S~'^ x)dx -^jT^ const ]^ A~'^{S^x)dx (4.3) 

j = -oo 

which is, clearly, a formal relation because pT tends to as T ^ oo. Note, however, that 
(4.3) implies that the ratios between prix) and pt{x'), with x^x' in the surface elements 
dx, dx' , are well defined, even in the limit as T — > oo. 

The equation (4.2) is for T large already a "good approximation" for the SRB distribution. 
It shows that the statistical averages should be computable by looking at a large part of Wq , 
namely at the finite approximation of the attractor A called above, and by imagining 
it coaied with a density pt{x), and then using (4.2). 

The existence of the limit as T ^ oo of (4.2) can be seen by remarking that the limit 
can in fact be written as an integral over phase space, in spite of the fact that pt{x) tends 
manifestly to as T — > oo. For, when T ^ oo what really matters is the amount of mass 
ending up inside a generic little square E in the phase space C, with center xe- Since E will 
be cut many times by 5^ A we can imagine that the various "pieces" of S'^ IS. intersecting 
E are piled up "vertically" in E: the figure below shows a picture for the simple case (4.1). 




Fig. 1: The parallel lines are intersections of the finite approximation for the attractor, S'^A, with the 
set E, represented by a square. The ^, r} axes are "parallel" to the unstable and stable manifolds Wq , 

respectively. Each of them is coated, eventually, by the image S'^U of U which gives them a thickness, (not 
shown). In the case of the map (4-1) tf^^ parallel lines are generated, if one moves on Wq away from O, 
in the following typical order: from bottom to top first one draws, successively, the lower line of each pair; 
then one draws the second, then one should draw a third series of lines above the second and keep going in 
this way until the endpoints of S"^ A are reached. For T — > oo the parallel lines fill densely E. 
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If iit{E) is the total mass initially in U ending up in E after time T, we can rewrite (4.2) 
as: 



Y,1^t{E)F{xe) 



(4.4) 



E 



provided E is so small that we can neglect the variation of F inside E, and that the _E's 
pave the phase space. Suppose we set up a coordinate system in the small box E (which 
is a box with full dimension 2D) so that the "horizontal" coordinates are called ^ and the 
"vertical" ones -q. A point in n A is denoted x{£,,ri) and the surfaces of constant 77 arc 
connected surface elements of the unstable foliation W" in E (a foliation of a set is a family 
of disjoint connected surfaces whose union is E), while those of constant ^ are connected 
surface elements of the stable foliation W. 

Then we see that (4.2), before the limit as T ^ 00 is taken, can be expressed as a sum 
over the connected parts of the surface^° S'^A that fall in E (the parallel lines in the above 
figure)). If {r]j} represents the j-th line then we can write (4.2) as: 



where c?^ denotes an area element on {rfj}. In the limit T 00, while PTi^^Vj) tends to 
the number of lines {rjj} tends to infinity and the sum over the surface elements of 5"^ A 
that cross E should converge to an integral over 77 and ^ with respect to some measure 
'p^{^)v{dri)d$_ with both the density 75''(^) along the unstable manifold and the measure v 
well defined. The measure v will give us the detailed information on how the various pieces 
{layers, or lines in fig. 1) of Wq intersecting E pile up and the distribution of the gaps 
between them in E,^^ hence on its fractal nature; on the other hand p''(C) will be a function 
such that the ratios J 'p^iCldE,/ J J/' (^)rfC should tell us the ratio of the masses of U ending 
up near the pieces of unstable manifold passing through r/ and 77', inside E, which should be 
well defined, by (4.3), in the limit T — > +00 as argued above. 

§5 The thermodynamic analogy. 

In this section we describe theoretical difiiculties with the heuristic analysis of §4 and with 
a mathematical proof of the existence of the limit (4.2). The solution to the difficulties that 
will be pointed out necessitates the introduction of more refined ideas and eventually the 
use of Markov partitions. We first point out the difficulty. 

If the analysis of the previous section is correct, i.e. if p''(0 really exists wc should 
be able to "calculate" it, at least formally. While it is evident that the function p''(^) is 

^ Which, wc recall, rcprcsnts the finite approximation to the attractor defined above. 

^ Note that "gap" here does not mean an actually empty region: since Wq is dense in E (by (C), §2) there 
can be no open regions in E which are not crossed by one connected part of Wg. In general, one should 
think of as supported by a dense "Cantor set" . 




/ p^{Oduu{drj)F{x{^,ri)) (4.5) 
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defined up to a constant for each rj such that (^,?7) G Wq, (because of formula (4.3) and 
the comment following it), it is much less evident that p''(^) behaves reasonably regularly 
in (^,r?) eEnW^. 

Clearly the right hand side of equation (4.3) can be used to compare the values of p''(0 
and o{'p^{^') if a; = (^, rj) and x' = (^', rj) are points of EdWq with the same r] (note that in 
such case only few of the factors in (4.3) have ratios really diferent from 1, because S~^x and 
S~^x' approach O exponentially fast and start close on a connected part of Wq). However 
if we compare p^i^) and (^) with the same ^ and {^,r]),{^,ri') G ECi Wq we run into the 
difficulty that the distance between (^, r)), (^, rj') m,6asured along Wq may be extremely long, 
in fact as long as we please by varying the two points on the surface ^ = constant. Therefore 
the function p^{^) might vary quite irregularly in E. Hence we see that the existence of a 
limit in (4.5) is not so obvious even though (4.3) provides immediately an expression for the 
ratios of the limit density on the set Wq D E. 

We m,ust find an alternative way to control the variations of such a function in _E n Wq in 
order to argue that p''(^) is well defined. The number of connected components of Wq in E is 
denumerable and we cannot expect that the SRB distribution is supported by a denumerable 
set of du dimensional surfaces. Hence wc arc in a position similar to when attempting to 
define the integral of a continuous function over a segment from the knowledge of the function 
at the rational points on the segment: this is possible only if the function is not too wildly 
changing from point to point. 

The resolution of this difficulty proceeds in two steps. First we push the analysis of the 
variability of p''(^) just given somewhat further to arrive at the equations (5.3), (5.4) below 
which are useful to illustrate the development of the thermodynamic analogy that gave rise 
to the thermodynamic formalism. This will enable us, in §6, to discuss the proper solution 
to the problem of the existence of the limit (4.5) and of the function p^(^), based on this 
thermodynamic analogy. 

Let X = {^,T]),x' = {£,,ri') be points of n Wq and let and d^' be two infinitesimal 
surface elements in £^ n Wq at different heights 77' and 77, corresponding to each other, in 
the sense that the stable manifolds through C W" intersect the unstable manifold W", 
exactly on d^', see fig. 2 below. Then the masses 'p^{^)d(_ on the segment d^ and on (O'^C 
have a ratio that can be computed by using (4.3) and by remarking that the ratio of the 
areas d^/d^' is: 



d^ _ \s-HsdO\ _ \s-^s^dO\ 

d^' \S-'^{Sd^')\ |5-2(5X')I 



\S-^{S^^dO\ 
\S-^{S^d^')\ 



(5.1) 



hence: 



Au'iS^C')) \S^d^'\ 



(n 



(5.2) 
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See Fig. 2 for an illustration: 




Fig. 2: The two lines at constant r] are at height rj and t]' respectively. The two infinitesimal segments 
and d^' correspond to each other as they are crossed by the same set of stable manifolds, the extreme two 
of which are drawn as dashed lines. 



But M^fx^ ^ because the two segments approach each other, while greatly and 

chaotically erring towards oo on Wq at the exponential speed of the expansion rates. 

Hence by combining (5.2) for M oo and (4.3), wc sec that the ratio of the measures 
P^iO^/P^ (^')<^^' corresponding intervals d^, d^' near the corresponding points x,x' ^ E 
with local coordinates (^, r?) and {^,r]'), is simply: 

This shows that the SRB distribution fi can be formally given by attributing to the "points" 
on the unstable manifold of A a weight given by: 

const H A-^{S-^x) = const e~^-'^^°^^"^^~'''^ (5.4) 

— oo 

or by a density on Wq given formally by the product in (4.3). Such statements should be 
interpreted in the same way in which one interprets statements like: "the one dimensional 
Ising model with nearest neighbour interaction attributes to the spin configuration {(Ti)??^ 
the probability" : 

const e S-oc ■^'^''^'+^ (5.5) 

The formal expression (5.4) must, therefore, be interpreted as a limiting statement. The 
function h{x) = logA„(a;) in (5.4) plays the same role as Jaoai in the Ising model in (5.5) 
and the appropriate way of understanding (5.4) is, as we have in fact discussed, as a limit 
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of (4.3). The important realization of Sinai, [S2], was that the equation (5.4), via (5.5), had 
a close analogy in statistical mechanics. "'^^ 

This remark led Sinai to his general theory of Markov partitions (see below) which are 
the main technical tool that is used to show mathematically that the limit of (4.2) or (more 
precisely) the limit in (4.5) really exist and for describing its general properties in satisfactory 
detail. 

%6 Coarse graining and Markov partitions. 

The above discussion has, as already said at the beginning of §5, just heuristic value as it 
has not led to a really usable formula for p''(^), but just to a few relations that such function 
must obey when evaluated on Wq. 

The solution lies in a stricter interpretation of the thermodynamic analogy (i.e. the 
similarity between (5.4) and (5.5)). To understand the rigorous solution (given in [S2]) 
to the problem of showing the existence of the limit (4.5) and the existence of and 
v{dvi), one has to introduce the concept of a "parallelogram" and of a Markov partition £ 
of the phase space C into parallelograms. This can be ultimately related to the problem of 
constructing a good division of the phase space in cells [i.e. a "good" coarse graining) so 
that the evolution can be correctly represented as a cell permutation, without "distorting" 
the hyperbolic nature of the motion (for such an interpretation of what follows see [Ga4]). 

A parallelogram will be a small set with a boundary consisting of pieces of the stable and 
unstable manifolds joined together as described below. The smallness has to be such that 
the parts of the manifolds involved look essentially "straight" : i. e. the sizes of the sides 
have to be small compared to the smallest radii of curvature of the manifolds W" and WJ, 
as X varies in C. 

Therefore let (5 be a length scale small compared to the minimal (among all x) curvature 
radii of the stable and unstable manifolds. Let W^'^,W^'^ be the connected parts of W^, 
containing x and contained in a sphere of radius 6. 

Let us first define a parallelogram E in the phase space C, to be denoted by A" x A^, with 
center x and axes A", A^ with A" and A*^ small connected surface elements on W!^ and 
containing x. Then E is defined as follows. Consider ^ £ A" and rj ^ and suppose that 
the intersection ^ x r] = W^'^ D W^'^ is a unique point (this will be so if 6 is small enough 
and if A", A® are small enough compared to S as we can assume, because the stable and 
unstable manifolds are "smooth" and transversal, see footnote ^). 

The set E = A" x A* of all the points generated in this way when ^, rj vary arbitrarily in 
A", A" will be called a parallelogram (or rectangle), if the boundaries 0A",9A* of A"and 
A" as subsets of and W.^, respectively, have zero surface area on the manifolds on which 
they lie. The sets duE = A" x dA'^ and dsE = OA"- x A** will be called the unstable or 
horizontal and stable or vertical sides of the parallelogram E. 



Hence the name of "thermodynamic formalism" given by Ruelle to the mathematical theory based on the 
above point of view, [Rl]. 

This is only approximately true because they are Holder continuous with some positive exponent, related 
to the gap between the positive or negative Lyapunov exponents and 0. 
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A*- 

X A' 



Fig. 3: The circles are a neighborhood of x of size very small compared to the curvature of the manifolds; the 

first picture shows the axes; the intermediate picture shows the X operation and W^'^ , W^'^ (the segments 
through rj and g have size 5); the third picture shows the rectangle E with the axes and the four marked 
points are the boundaries 9 A" and dA". The picture refers to a two dimensional case and the stable and 
unstable manifolds are drawn as flat, i.e. the A's is very small compared to the curvature of the manifolds. 



Consider now a partition £ — [Ei, . . . , Ejsf) of C into A/" rectangles Ej with pairwise 
disjoint interiors. We call 9„£ = LijduEj and ds£ = UjdgEj-. these are called respectively 
the unstable boundary of £ and the stable boundary of £, or also the horizontal and vertical 
boundaries of £, respectively. 

We say that f is a Markov partition if the transformation S acting on the stable boundary 
of £ maps it into itself (in formula this is: Sds£ C ds£) and if likewise the map acting 
on the unstable boundary maps it into itself {S~^du£ C du£)- 

The actual construction of the SRB distribution then proceeds from the important result 
of the theory of Anosov systems expressed by what we shall call "Sinai's first theorem" : 

Theorem: every transitive Anosov system admits a Markov partition £, [S2]. 

The above theorem is the first step towards a controlled version of the heuristic arguments 
given above and towards a usable form of equation (4.5) based on a suitable interpretation 
of (5.3). It can be extended to imply the existence of more special Markov partitions: for 
instance to show the existence of Markov partitions with any one of the following three 
properties (the last shows that the first two can be realized simultaneouly and will play a 
key role in our analysis): 

1) The construction of £ can be done, [Ga3], so that the horizontal axes of Ej all lie on 
Wq (and the vertical on Wq) and their union is a set that can be obtained from a single 
small connected surface element A of Wq (resp. A of Wq) containing O by dilating it 
with a high iterate S"^ of the time evolution S. In other words the union UjA" of the 
horizontal axes of the parallelograms Ej e £ can be regarded as a good finite approximation 
to our attractor A, because it has the form S^A. with A a connected surface element of the 
unstable manifold Wq, containing O. Likewise the union of the stable axes can be regarded 
as a large connected part of the stable manifold Wq. 

2) If the reversibility property holds it is clear that i£ is also a a Markov partition. This 
follows from the definition of Markov partition and from the fact that reversibility implies: 

= iWl (6.1) 
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The definition of a Markov partition also implies that the intersection of two Markov par- 
titions is a Markov partition, hence it is clear that there are Markov partitions £ that are 
reversible in the sense that £ = i£. 
3) Furthermore one can construct a Markov partition £ which is reversible and at the 

same time verifies the property 1) above, [Ga3]. 
Here we shall use Markov partitions that verify property 3) above. 

In order to formulate Sinai's second theorem, which gives an expression for a controlled 
approximation to the SRB distribution, we consider the partition £t = r\^j,S^^£ obtained 
by intersecting the images under , j = —T, . . .,T o{£. Then £t is still a Markov partition 
and it is time reversal invariant if £ is. We now construct a probability distribution that 
we regard concentrated on the finite approximation At to the attractor, consisting of the 
union of the horizontal axes of £t, (see remark 1),2) above), and equal, with the notations 
of rem ark 1 above, to 8"^+^ A = At- 

We can visualize the small parallelograms forming £t as a lattice of parallelograms: two 
parallelograms adjacent and on the same vertical strip will have horizontal axes that cor- 
respond to each other in the same sense that the close horizontal surface elements used in 
deriving (5.3) correspond. Therefore if wc attribute to the horizontal axis of the parallelo- 
gram Ej in £t with center Xj a weight equal to Y\^fl2-T/2^u^^^^''^i) ^^'^^ ratios 
of the weights of corresponding surface elements automatically realize an approximation of 
the product in (5.3). We take, of course, t ^ T so that the size of each parallelogram is 
so small that the weight we attribute to each does not depend on which point of Ej we 
regard as a center and that no essential ambiguity arises as to which weight to attibute 
to a parallelogram.^^ Note that the above weight A„.^(a;) is the inverse of the expansion 
coefficient of the map S"^ as a map of W^_^/2^ to W^^/^^ (between S~'^/'^x and S'^/'^x), i.e. 

r/2-l 

n ^u{S'^) (6.2) 

j=-r/2 

A similar quantity can be defined by regarding as a map of Wg_^/2^ to WJ^/a^. 

The construction thus generates a probability distribution which, by the above analysis, 
verifies (5.3) more and more exactly as r — » oo. Hence this analysis suggests the following 
theorem (which we shall call "Sinai's second theorem", [S2]): 

Theorem: If (C, S) is a transitive Anosov system the SRB distribution fi exists and the fi 
average of a smooth function F is: 

I n{dx)F{x)= lim hm [nTr{dx)F(x) (6.3) 

where, with the above notations, xj is a point in Ej G £t . 

The above /ZT,r as defined by the middle ratio in (6.3), can be taken as a "concrete" 
procedure to follow in approximating the SRB distribution. 



The size of the pajrallelograms of £t is cleajrly decreasing as e , at least, if A is the spectral gap, see 
footnote ^ above. 
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In the case of equilibrium under assumption (C) in §2 the distribution fj, in equation (6.3) 
can be shown to coincide with the micro canonical ensemble, as already mentioned, [Sl].^^ 



§ 7 Application. 

The chaotic hypothesis can be taken as an extension of the ergodic hypothesis for equi- 
librium statistical mechanics to systems in nonequilibrium stationary states (conservative 
or dissipative). In the equilibrium {i.e. conservative) case it implies the ergodic hypothesis 
{but it is stronger) and hence the microcanonical distribution, which we know how to use in 
order to draw physical consequences. 

It is therefore legitimate to ask whether the chaotic hypothesis and the ensuing SRB 
distribution have any predictive value of their own. Just as the ergodic hypothesis implies 
the well tested classical thermodynamics, the new hypothesis should imply, for example, 
irreversible thermodynamics of nonequilibrium stationary states, without the necessity of 
solving the equations of motion. It is not clear that this is so. 

However there are already some experimental results that offer support to the chaotic 
hypothesis, since one can understand their outcome by using it. 

Here we examine, in particular, one experimental result, [ECM2], which the authors al- 
ready attempted to explain by relating it to our chaotic hypothesis. Some of the data in 
[ECM2] require, to be unambigously understood, the discussion in [ES]. We shall take the 
viewpoint of the preceding sections to make more precise, and detailed, the argument in 
[ECM2], modifying it to some extent in order to put it on a more mathematical basis. 

In the context of this paper the experiment [ECM2] deals with model 2 in §3 and measures 
the entropy production rate {i.e. phase space contraction rate, see §3, (3.5)) as seen on a 
stretch of time r, short compared to the duration of the experiment T, and repeating the 
measurement T/r times. We emphasize that this is an experiment on a system far from 
equilibrium. 

Calling D(Tr{x) the entropy production rate measured on the motion originating at S^'^^'^x 
and observed r units of time (we take r even for simplicity), we define it, see (3. 4), (3. 5), by: 



where {o)_^ is the average in the future of a{S^x), which is a constant almost everywhere in 

^ One should not be disturbed by the fact that this is a rigorous mathematical theorem only for Anosov 
systems or, more generally, for "axiom A" attractors, [Rl]: one should not forget that the microcanonical 
ensemble is also lacking a mathematical justification in equilibrium theory. In fact the only equilibrium case 
in which one can prove the ergodic hypothesis is for the Lorentz gas (i.e. the billiards), [SI], [BSC], with 
N = 1: in such a case (N = 1, = and ip" a suitable hard core potential, i.e. a triangular lattice of hard 
disks) the present point of view can also be shown to hold in the presence of dissipation, [CELS]. 
^ The reader should not mind that the symbol for the integer T is sometimes also used for the absolute 
temperature. 




t/2-1 



(7.1) 



J=-r/2 
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phase space with respect to /zo-random choices of initial data.^'' 

The total entropy production while the phase space point x evolves between S^^/'^x and 
S'^l'^x is obtained by multiplying (7.1) with the time elapsed during r such timing collisions. 
For simplicity we think that the time interval to between the timing collisions is constant. 
Note that x is the middle point of the segment of a trajectory of (discrete) time length r, 
defining the fluctuation ar{x) in (7.1). 

It is perhaps important to stress that (cr)_,_ is very different from the limit as t ^ +00 of 
arix) in (7.1): in fact the latter (by the time reversal symmetry) vanishes, while the former 
is positive as follows from numerical evidence, see §5, and as assum,ed in (A) in ^2. 

The experiment divided the ai--axis into small intervals /o,/±i,... and measured the 
quantity ar{x) building a histogram counting how many times the fell into the interval 
Ip (where a,-(x) = p). Obviously we expect a distribution 7ri-(p) centered around an average 
which is, (see (3.4) and (7.1)), exactly 1. The result for TTrip) can be found in fig. 1 of 
[ECM2], for one (rather large) value of r and one of 7 and N = 56. 

A second experimental result is for 11^ (p) = — 2jvrto(o-) t^^(-p) ' essentially for the 
logarithm of the ratio of the probability that a,-(a;) = p to that of a,- (a;) = —p. The result, 
fig. 2 of [ECM2] is, for the rather large value of t considered, a remarkably precise straight 
line for 11^ (p) as a function of p, i.e. n^(p) is a linear function of p. 

A third experiment shows that the slope of this line as a function of r verifies, even for 
large deviations, the relation of proportionality to r for r large (fig. 3 of [ECM2]). 

The results are rather precise with apparently little margin for errors, hence one has to 
find a theoretical reason that the probability distribution of D(Jt{x) has the form: 

7r,(p) dp'^=P{ar e{p,P + dp)) = e-^^^P^+^^Pdp 

7rr(-p) 

for a suitably chosen constant C and a suitable even function C(p) with minimum at |? = 1 and 
with the argument of the exponential correct up to, apparently, p, r independent corrections 
(see fig. 3 of [ECM2]). In [ECM2] a theoretical argument is presented which leads to 
2C = Dt{){(T) j_, if to is the average time between timing events. 

We are now going to show, and this is our main technical result (and a theorem under 
assumptions (A,B,C) of §2), what we call a fluctuation theorem: 

Fluctuation theorem: Let {C,S) verify the properties (A,B,C) of%2, (dissipativity, reversibil- 
ity and chaoticity). Then the probability iTt{p) that the total entropy production DTtoCTrix), 
(7.1), over a time interval t = rto (with to equal to the average time between timing events) 
has a value Dt{a)j^p verifies the large deviation relation: 

= e^*<'^>+^' (7.3) 

t^t{-p) 

This because the SRB distribution verifies the extended zeroth law, (2.1), which says that the averages 
are, with /io-probabiUty 1, independent of the initial data. 



or 

(7.2) 
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with an error in the argument of the exponential which can be estimated to be p, r-indepen- 
dent. 

This means that if one plots the logarithm of the left hand side of (7.3) as a fimction of p 
one observes a straight line with more and more precision as r becomes large (in agreement 
with figure 3 in [ECM2]). 

Remark: since the above theorem is deduced under the assumptions (A,B,C) only, the 
result (7.2) will apply as well to the models 1,3,4,5. This gives a parameterless prediction of 
the outcome of several numerical experiments similar to the one described above. 

The main ideas for the proof, [CG], [Ga3], [Ga2], of the above theorem are the following. 

The probability that a^(a;) G Ip over the probability that arix) £ /_p is, if one uses 
the notations and the approximation /Ut,t to /x described at the end of §4 (see (6.3)) with 
F{x) = ar{x): 

T,j,a^{xj)=p K,r{Xj) ^^^^ 

where Au,t{x) is the jacobian determinant of 5"^ as a map of Wq into itself, evaluated at 
the point e {i.e. as a map between S-'^/'^x and S^/^x, see (6.2), (7.1)). 

Since /iT.r in (6-3) is only an approximation at fixed T, r an error is involved in using 
(7.4)- It can be shown that this error can be estimated to affect the result only by a factor 
bounded above and below uniformly in r, p, [CG] , [Ga3] , [Ga2] . This is a remark technically 
based on the thermodynamic analogy pointed out in (5.4), (5.5). 

We now try to establish a one to one correspondence between the addends in the numerator 
of (7.4) and the ones in the denominator, aiming at showing that corresponding addends 
have a constant ratio which will, therefore, be the value of the ratio in (7.4). 

This is possible because of the reversibility property (B), §2. Let x € A, then ix e iA. By 
using the identity S~'^{S'^x) = x, the identity S~'^{iS~'^x) = ix (time reversal) and (6.1), 
we deduce the relations:^* 

ar{x) = -ar{ix), A„,^(ia;) = A^l{x) (7.5) 

which are identities, see [ECM2],[CG],[Ga3]. The first equality in (7.5) is obvious as in all 
the cases considered the i operation changes the sign to a{x), the rate of change of the 
phase space volume. The second equality in (7.5) is also easy to check: in fact let /3 be a 
surface element on W" aroimd S^'^^'^x and let /3' = S'^f3 be its S'^ image around S'^^^x: 
then Au,t)x) = Applying i to (3 and /?' one obtains surface elements i/? and i/3' on 
with the same area as /3 and /?' (because i is an isometry) around respectively S'^^^ix and 
S~'^/'^ix so that the expansion rate As_^(ia;) is = = A„_^(a;). 



The key remark is that time reversal i maps Ej into iEj and at the same time changes the horizontal 
surface elements of Ej into the vertical ones of iEj and the vertical surface elements of Ej into the horizontal 
of iEj, see (6.1). 
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The ratio (7.4) can therefore be written simply as: 



Eb, ,a, (x,-) = -p K!t {Xj ) ^E, ,a. {x,)=p Kr i^Xj ) 

where Xj G Ej is a point in E'j . In deducing the second relation we make us of the existence 
of the time reversal symmetry i and of (7.5). 
It follows then that the ratios between corresponding terms in the ratio (7.6) is equal 

to A„ ^(a;)Ag ^(x). This differs from the reciprocal of the total variation of phase space 
volume over the time r between the point S^'^^^x and S^^^x only because it does not 
take into account the ratio of the sines of the angles 'd{S"^^'^x) and ■d{S'^^^x) formed by 
the stable and unstable manifolds at the points S~'^/^x and S'^/'^x, see footnote ^. But 
A„ ^(a;)Ag ^(a;) will differ from the actual phase space contraction mider the action of S*^, 
as a map between S^'^/'^x and S'^^^x, by a factor that can be bounded between and B 
with B = maxa;^j;/ | which is finite by the transversality of the stable and unstable 

manifolds. 

Now for all the points Xj in (7.6), the reciprocal of the total phase space volume contraction 
over a time to is e~°^^'^^'-'^'^^+*°^''^, which (by the constraint imposed on the summation labels 
= p) equals e-^*'''"<'">+^'. Hence the ratio (7.4) will be e^*'""<''>+^', to leading order as 
N,T ^ cxD, proving (7.3), with 2C = D{a)^tQ. It is important to note that there arc two 
errors ignored here, as pointed out in the discussion above. They imply that the argument 
of the exponential is correct up to p, t independent corrections (which are in fact observed 
in the experiment as fig. 3 of [ECM2] shows). 

The p independence of the coefficient of C in (7.2) is therefore a key test of the theory 
(and it should hold with corrections of order 0{t~^)). 



§5 Outlook. 
We end with a number of remarks. 

1) The interest of our discussion in §7 is not, of course, the fluctuation theorem which is 
essentially proved there (for a formal proof see [Ga2] and [Ga3] ) , but in the clarification of the 
meaning of the properties (A,B,C) mentioned in §2. Furthermore it is interesting that our 
chaotic hypothesis, §2 does have some concrete and experimentally verifiable consequences 
(verified here in the case of model 2): such consequences have the remarkable feature of 
being predictions without free parameters, hinting that the hypothesis might have a quite 
general validity. One cannot be too demanding on the matter of mathematical rigor: one 
should not forget that the ergodic hypothesis is far from being proved either, particularly in 
the generality one would want. 

2) The fluctuation formula (7.3) holds also for the models 1,3,4,5 because the fluctuation 
^ In the previous paragraph and in the paragraph following (7.4). 
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theorem applies to all such models (see remark following the theorem): but the numerical 
experiments do not seem to exist yet. 

3) The pairing property (D) and the smooth distribution of the Lyapunov exponents (E) 
have been used here only to get some intuition and to visualize the hyperbolic nature of the 
attractor and the equality of the dimensions of the stable and unstable manifolds. It seems 
interesting to perform numerical experiments to try to investigate better, at least in the 
systems that we are considering, if the density function foo{x) is really positive at a; = 0, 
see (E), §2, as the numerical results seem to suggest in some cases, [LPR], [ECMl], [SEM]. 

4) Note that the fluctuation theorem ((7.3)) applied to model 5 leads to an interesting 
consequence on the large deviations properties of the magnitude of the energy dissipation s 
in turbulent flows. In this case we take, for simplicity, the kinetic energy DykT{x) of the 
viscous modes to be a non fluctuating quantity equal to DykT. Then the random variable 
p associated with a{x) = ^ ^^^^^^^ in the fluctuation theorem of §7 is just proportional to 
the average over a time interval t of the energy dissipation rate e. This is a variable that is 
assumed to be constant in the Kolmogorov-Obuchov theory: what we say here is that it is 
in fact a fluctuating quantity and we predict (on the basis of the fluctuation theorem of §7) 
that the time average (e)^, over a time interval t, of s{x),i.e. {s)^ = tDy{e)_^p is such that 
its probability distribution ttj {p) verifies the linear large deviation law. 

n^^^o,[.M/M-P)]^i (8.1) 

up to corrections of order 0{t^^). If T = T{x) has to be regarded as a fluctuating variable 
then (8.1) must be regarded as a property of the fluctuations of the entropy production 
rate a{x) = rather than of the energy dissipation, (with some obvious modifications, 

e.g. {e)JkT 

5) Concerning the particularity of the gaussian thermostat we think, see [Gal],[C], that 
there should be, also in nonequilibrium, several equivalent ways of describing the same 
stationary distribution corresponding to different jj. and to different physical ways of reaching 
the stationary state, at least in the thermodynamic limit. And it may well be that the 
gaussian thermostat turns out to be equivalent to other models of thermostats, which could 
be described by rather different attractors. For instance a stochastic thermostat, in which 
a particle colliding with the wall is scattered with a maxwellian distribution at a given 
temperature, will certainly be described by a statistics /x which is absolutely continuous 
with respect to the Liouville distribution.^" In the thermodynamic limit this might just give 
the same result as obtained with a statistics which, for finite N, is on a fractal attractor. This 
mechanism is like the one realized by the microcanonical and the canonical ensembles (the 
first is concentrated on a set of configurations which has zero probability with respect to the 
second, as long as N < oo). This is clearly a question that requires further investigations. 

6) One can also regard the gaussian thermostat as a device to eliminate some trivial Lya- 

Note that a stochastic model of thermostat is described by a stochastic differential equation and therefore 
our discussion does not apply without some major modification. 
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punov exponents. For instance in model 4 we could simply not introduce the gaussian 
constraint that the total energy is constant (which led to the ctoP^ terms): we believe that 
physically the system would then still behave in the same way, for large N, L. But without 
such a constraint we could not assume the phase space to have 4A'' — 3 dimensions, because 
H would not be rigorously constant. Wc expect, however, that in such case the energy H is 
approximately constant, in fact more and more so as L,iV ^ oo with NL~'^ = n constant. 
Thus the variability of H probably leads to a zero Lyapunov exponent making the chaoticity 
assiimption manifestly invalid. But this would be so only in a somewhat trivial way as its 
violation is due to a zero Lyapunov exponent associated with a variable which is "almost" a 
constant of the motion. Hence it is natural to fix the value of the energy rigorously a priori 
by a constraint (realized via a minimum constraint principle, like Gauss' principle) and so 
dispose of the extra (or very close to 0) Lyapunov exponent, recovering then again a situ- 
ation in which the system is strictly chaotic. Such a point of view can be extended to cover 
cases in which hyperbolicity is not valid because of the existence of quasi exact conservation 
rules. An example is in fact model 2 in which the variables yj,Pyj can be replaced by yj,pyj 
thus turning P^jX^ into variable quantities: their variability is however clearly due to the 
special (vertical) boundary conditons used and it should therefore not matter whether they 
are kept rigorously constant or not, in the limit of TV, L ^ oo. The dynamics can be modified 
by turning such quantities into exact conservation laws and the new dynamics should be 
indistinguishable from the previous one in the thermodynamic limit. Another example is 
provided by the constraints imposed on model 3 to achieve that the horizontal momentum 
is conserved. 

7) Like for the ergodic hypothesis in equlibrium, the range of validity of the chaotic hypothe- 
sis for nonequilibrium stationary states is not known; the more complicated the nature of the 
latter states, maintained in the presence of external fields or special boundary conditions, 
makes this case even more difficult. 
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